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Abstract 

New experimental results on polarized structure functions, cross sections 
for e^p neutral and charge current reactions and v iv) charge current on 
isoscalar targets are compared with predictions using the statistical parton 
distributions, which were previously determined. New data on cross sections 
for Drell-Yan processes, single jet in pp collisions and inclusive 7r° produc- 
tion in pp collisions are also compared with predictions from this theoretical 
approach. The good agreement which we find with all these tests against 
experiment, strenghtens our opinion on the relevance of the role of quantum 
statistics for parton distributions. We will also discuss the prospects of this 
physical framework. 
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1 Introduction 



Deep-inelastic scattering (DIS) of leptons and nucleons is, so far, our main 
source of information to study the internal nucleon structure, in terms of 
parton distributions. Three years ago we proposed to construct, in a 
unique way, the unpolarized and the polarized parton distributions, using a 
simple procedure, inspired by a quantum statistical picture of the nucleon, 
in terms of Fermi-Dirac and Bose-Einstein functions. An important feature 
of this new approach lies into the fact that the chiral properties of pertur- 
bative quantum chromodynamics (QCD), lead to strong relations between 
quark and antiquark distributions. As a consequence the determination of 
the best known unpolarized light quarks (u, d) distributions and their cor- 
responding polarized ones (Au,Ad), allows to predict the light antiquarks 
distributions (u, d, Au, Ad). Therefore our approach has a strong predictive 
power, in particular, the flavor-asymmetric light sea, i.e. d > u, which can 
be understood in terms of the Pauli exclusion principle, based on the fact 
that the proton contains two u quarks and only one d quark It is also 
natural to anticipate that the signs of Au and Ad are the same as Au and 
Ad, respectively. One more relevant point to recall, is that all these parton 
distributions were determined in terms of only eight free parameters, which 
is indeed remarkable. 

More recently we compared j3] our predictions with some new unpolarized 
and polarized DIS measurements obtained at DESY, SLAC and Jefferson 
Lab. and they turned out to be rather satisfactory. Therefore, in order to 
strengthen the relevance of this physical picture, we carry on the comparison 
with data from a much broader set of processes, including new DIS results 
and also hadronic cross sections. 

The paper is organized as follows. In the next section, we review the 
main points of our approach for the construction of the statistical parton 
distributions and we recall their explicit expressions. In section 3, we discuss 
in more details the predictive power of our approach in connection with 
some simple mathematical properties of the Fermi-Dirac expressions and the 
numerical values we found for the free parameters. It allows us to clarify 
the x behavior of the quark distributions, where it is known from the data, 
but also to foresee some specific behaviors, in so far unexplored regions, 
for example in the high x domain. In section 4, we consider e ± p neutral and 
charged current reactions, whereas section 5 deals with u(u)p charged current 
reactions. Section 6 concerns Drell-Yan processes, while section 7 deals with 
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inclusive single-jet production in pp collisions and inclusive 7r° production in 
pp collisions. We give our final remarks and conclusions in the last section. 



2 The quantum statistical parton distribu- 
tions 

The light quarks q = u, d of helicity h = ±, at the input energy scale 
Qq = 4GeV 2 , are given by the sum of two terms |T], a quasi Fermi-Dirac 
function and a helicity independent diffractive contribution, common to all 
light quarks 

h( 2N _ A K xb 

xq (x, Q ) _ x ^_^ + i + exp(x/a;) + 1 - [1) 

Here Xq q is a constant, which plays the role of the thermodynamical potential 
of the quark q h and x is the universal temperature, which is the same for 
all partons. We recall that from the chiral structure of QCD, we have two 
important properties, allowing to relate quark and antiquark distributions 
and to restrict the gluon distribution |U E] : 

- The potential of a quark q h of helicity h is opposite to the potential of 
the corresponding antiquark q~ h of helicity -h 

K = -X$ . (2) 

- The potential of the gluon G is zero 

X 0G = . (3) 

Therefore similarly to Eq. (|TJ). we have for the light antiquarks 

h ( A(X - fe )-V fc Ax 1 

xq {x, Q ) exp[(a; + x _ h)/s] + i + exp(x/ _ ) + 1 . (A) 

Here we take 2b for the power of x and not b as for quarks, an assumption 
which was discussed and partly justified in Ref. [lj. 

Concerning the unpolarized gluon distribution, we use a quasi Bose-Einstein 
function, with zero potential 

xG(x,Ql) = — - . (5) 

exp(x/x) — 1 
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This choice is consistent with the idea that hadrons, in the DIS regime, are 
black body cavities for the color fields. It is also reasonable to assume that for 
very small x, xG(x, Ql) has the same behavior as the diflractive contribution 
of the quark and antiquark distributions in Eqs. (JTJ) and (jlj), so we will take 
be = 1 + b. We also need to specify the polarized gluon distribution and we 
take 

xAG(x,Q 2 ) = , (6) 

assuming a zero polarized gluon distribution at the input energy scale Ql . 
For the strange quarks and antiquarks, s and s, given our poor knowledge 
on their unpolarized and polarized distributions, we take 2 

1 

xs(x, Ql) = xs(x, Ql) = ~[xu(x, Ql) + xd(x, Ql)} , (7) 

and 

xAs{x, Ql) = xAs(x, Ql) = \[xAd{x, Q 2 ) - xAu{x, Ql)} . (8) 

This particular choice gives rise to a large negative As(x, Ql). Both unpolar- 
ized and polarized distributions for the heavy quarks c, b, t, are set to zero 
at Ql = 4GeV 2 . 

With the above assumptions, we note that the heavy quarks do not in- 
troduce any free parameters, likewise the gluons, since the normalization 
constant Aq in Eq. (jSJ) is determined from the momentum sum rule. Among 
the parameters introduced so far in Eqs. (0) and A and A are fixed by 
the two conditions u — u = 2, d — d = 1. Clearly these valence quark con- 
ditions are independent of b and A, since the diflractive contribution cancels 
out. Therefore the light quarks require only eight free parameters, the four 
potentials Xq u , Xq u , XQ d , X^ d , one universal temperature x, b, b and A. 

From well established features of the u and d quark distributions extracted 
from DIS data, we anticipate some simple relations between the potentials: 

- u(x) dominates over d(x), therefore one expects Xq u + X 0u > X^ d + X^ d 

- Au(x) > 0, therefore X^ u > X ~ 

- Ad(x) < 0, therefore X~ d > X+ d . 

So Xq u should be the largest thermodynamical potential and X^ d the 
smallest one. In fact, as we will see below, we have the following ordering 

Xl> X m ~ X 0u > X+ . (9) 

2 A strangeness asymmetry, s(x) ^ s(x), can be also obtained in the statistical approach 

CI 
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This ordering leads immediately to some important consequences for quarks 
and antiquarks. 

First, the fact that X^ d ~ X^ u , indicated in Eq. Q, leads to 

u-(x,Q 2 )<d-(x,Ql) , (10) 

which implies from our procedure to construct antiquark from quark distri- 
butions, 

u + (x,Q 2 Q )>d + (x,Ql) . (11) 

These two important approximate relations were already obtained in Ref . P , 
by observing in the data, the similarity in shape of the isovector structure 
functions 2xg ( f~ n \x) and F!f~ n \x), at the initial energy scale, as illustrated 
in Fig. n 3 - For 2xg ( f~ n \x) the black circles are obtained by combining 
SLAC ^T] and JLab ^2] data. The white circles, which extend down to 
the very low x region, include the recent deuteron data from COMPASS 13 
combined with the proton data from SMC [TU], at the measured Q 2 values 
of these two experiments 4 . The agreement with the curve of the statistical 
model is improved in this later case. The + helicity components disappear 
in the difference Ixgf Since this difference is mainly non- 

zero for 0.01 < x < 0.3, it is due to the contributions of u~~ and d~ (see 
Ref. P). 

Second, the ordering in Eq. Q implies the following properties for anti- 
quarks, namely: 

i) d(x) > u(x), the flavor symmetry breaking which also follows from the 
Pauli exclusion principle, as recalled above. This was already confirmed by 
the violation of the Gottfried sum rule [HI E] . 

ii) Au(x) > and Ad(x) < 0, which have not been established yet, given 
the lack of precision of the polarized semi-inclusive DIS data, as we will see 
below. One expects an accurate determination of these distributions from the 
measurement of helicity asymmetries for weak boson production in polarized 
pp collisions at RHIC-BNL (Hj, which will allow this flavor separation. 

By performing a next-to-leading order QCD evolution of these parton 
distributions, we were able to obtain in Ref. jT], a good description of a 
large set of very precise data on the following unpolarized and polarized DIS 

3 Notice that it differs from Fig. 1 in Ref. pQ, where we put incorrect scales, both on 
the vertical and the horizontal axes 

4 We have not included some corrections due to difference of the beam energies of 
COMPASS and SMC 
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Figure 1: The isovector structure functions Ixgf (x) (solid line from our 

statistical parton distributions) and F 2 (x) (dashed line from our statisti- 
cal distributions). Data are from NMC 0, SMC [TOJ, SLAC jTT], JLab [12] 
and COMPASS US]. 

structure functions F%' ' n (x, Q 2 ), xFg N (x, Q 2 ) and g\ ,d,n {x,Q 2 ), in a broad 
range of x and Q 2 , in correspondance with the ei<?/i£ free parameters : 

X + = 0.46128, X ~ = 0.29766, X~ = 0.30174, X+ = 0.22775 , (12) 

x = 0.09907, b = 0.40962, b = -0.25347, A = 0.08318 , (13) 
and three additional parameters, which are fixed by normalization conditions 

A = 1.74938, A = 1.90801, A G = 14.27535 , (14) 
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as explained above. Note that the numerical values of the four potentials 
are in agreement with the ordering in Eq. (jUJ), as expected, and all the free 
parameters in Eqs. (fT2l [TH|) have been determined rather precisely, with an 
error of the order of one percent. 

3 The predicting power of the statistical par- 
ton distributions 

We now try to relate the x dependence of the quark (antiquark) distributions 
to their specific expressions given in Eqs. (JTJ) and (JU) and to study the role 
of the different free parameters involved, according to their numerical values 
obtained in Ref. £Q. First, it is useful to note that, given the small value 
of A (see Eq. (fT3|) ). the diffractive contribution is less than 1CT 2 or so, for 
x > 0.1, but it dominates in the very low x region, when x « x, since 
b < 0. Therefore the strong change of slope of xu(x) and xd(x) at high x (at 
the input scale Ql and above), is related to the values of the corresponding 
potentials and is larger for u than for d, because of the ordering in Eq. Q. 
This is indeed what we observe in Fig. El at some rather high Q 2 values. This 
feature is not spoilt by the Q 2 evolution, which is also well described by the 
statistical quark distributions as shown in Fig. El where we compare with HI 
data. Another interesting point concerns the behavior of the ratio d(x)/u(x), 
which depends on the mathematical properties of the ratio of two Fermi-Dirac 
factors, outside the region dominated by the diffractive contribution. So for 
x > 0.1, this ratio is expected to decrease faster for X^ d — x < x < Xq u + x 
and then above, for x > 0.6 it flattens out. This change of slope is clearly 
visible in Fig. EJ with a very little Q 2 dependence. Note that our prediction 
for the large x behavior, differs from most of the current literature, namely 
d(x)/u(x) — > for x — > 1, but we find d(x)/u(x) — > 0.16 near the value 1/5, 
a prediction originally formulated in Ref. jTHj. This is a very challenging 
question, since the very high x region remains poorly known, as shown in 
Fig. |21 and Fig. El The typical behaviour of the Fermi-Dirac functions, falling 
down exponentially above the thermodynamical potential, which shows up 
in Fig. Q complies well with the fast change in the slope of gf{x) at high x, 
as shown in Fig. El 

Analogous considerations can be made for the corresponding helicity dis- 
tributions, whose best determinations are shown in Fig. El By using a sim- 
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ilar argument as above, the ratio Au(x)/u(x) is predicted to have a rather 
fast increase in the x range (Xq u — x, Xq u + x) and a smoother behaviour 
above, while Ad(x)/d(x), which is negative, has a fast decrease in the x 
range (X^ d — x, X^ d + x) and a smooth one above. This is exactly the trends 
displayed in Fig. El and our predictions are in perfect agreement with the 
accurate high x data 5 . We note the behavior near x — 1, another typical 
property of the statistical approach, also at variance with predictions of the 
current literature. The fact that Au(x) is more concentrated in the higher x 
region than Ad(x), accounts for the change of sign of (x), which becomes 
positive for x > 0.5, as first observed at Jefferson Lab [T2"] . 
For the light antiquark distributions (see Eq. (0J)), it is clear that in the very 
low x region (x < 10~ 3 ) the ratio d(x)/u(x) is ~ 1, since the diflractive con- 
tribution dominates 6 and it is an increasing function of x because the non 
diflractive term is larger for d(x) than for u(x). This natural expectation, 
d(x) > u(x) from the statistical approach, was already mentioned above and 
has been also confirmed by the E866/NuSea Drell-Yan dilepton experiment 
[20J, up to x = 0.15. For larger x, although the errors are large, the data 
seem to drop off in disagreement with our predictions (see Fig. 16 in Ref. PQ). 
This important point deserves further attention and we will come back to it 
in Section 6, when we will discuss Drell-Yan dilepton cross sections. This 
is another challenging point, which needs to be clarified, for example with 
future measurements by the approved FNAL E906 experiment [2J, to higher 
x values. 

We now turn to the antiquark helicity distributions. Since we predict Au(x) > 
and Ad(x) < 0, the contribution of the antiquarks to the Bjorken sum 
rule (BSR) [22| i s m our case 

0.022, at Q 2 = 5GeV 2 , which is not negligi- 
ble. The statistical model gives for the BSR the value 0.176, in excellent 
agreement with the QCD prediction 0.182 ± 0.005 and with the world data 
0.176 ±0.003 ±0.07 jTJ. It is also interesting to remark that Eq. (fTT|) implies 

Au(x) - Ad(x) ~ d(x) - u(x) > , (15) 

so the flavor asymmetry of the light antiquark distributions is almost the 
same for the corresponding helicity distributions. Similarly, Eq. (jlOj) implies 

It is worth mentioning that the Jefferson Lab points for the d quark are those of 
Ref. JSIi which have been moved down compared to those of Ref. 12 
6 Obviously, this is also the case for the ratio d(x)/u(x) 
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Figure 2: Statistical quark distributions xu(x, Q 2 ), xd(x,Q 2 ) as a function 
of x for Q 2 = 3000, 8000GeV 2 . Data from HI [B]. 



Au(x) - Ad(x) ~ u(x) - d(x) > . (16) 

By combining Eqs. (JT5|l and (|16p. we find a very simple approximate result for 
the BSR, namely ~ 1/6, a value compatible with the numbers quoted above. 
We also compare in Fig. our predictions with an attempt from Hermes 
to isolate the different quark and antiquark helicity distributions. The poor 
quality of the data does not allow to conclude on the signs of Au(x) and 
Ad(x), which will have to wait for a higher precision experiment. 
Finally we are coming back to the polarized gluon distribution AG(x,Q 2 ), 
which was assumed to be zero at the input scale Ql = 4GeV 2 (see Eq. (JHJ). 
It is interesting to note that after evolution, it becomes negative for Q 2 < Qq 
and positive for Q 2 > Qq. The results are displayed in Fig. |H]and are waiting 
for an improved experimental determination of this important distribution. 
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Figure 3: Statistical quark distributions c • xu(x, Q 2 ), c ■ xd(x,Q 2 ) as a 
function of Q 2 for fixed x bins. Data from HI [To] . 
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Figure 4: The quarks ratio d/u as function of x for Q 2 = 4GeV 2 (solid line) 
and Q 2 = lOOGeV 2 (dashed-dotted line). 
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Figure 5: Ratios (Au + Au)/(u + u) and (Ad + Ad)/(d + d) as a function 
of x. Data from Hermes for Q 2 = 2.5GeV 2 and a JLab experiment [T??] . 
The curves are predictions from the statistical approach. 
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Figure 6: g{(x, Q 2 ) as a function of x at fixed Q 2 = 3GeV 2 from the statistical 
approach. Experimental data from SLAC E143 [TP] . 
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Figure 7: Quarks and antiquarks polarized parton distributions as a function 
of x for Q 2 = 2.5GeV 2 . Data from Hermes [23]. The curves are predictions 
from the statistical approach. 
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Figure 8: The ratio AG(x)/G(x) as a function of x, for Q 2 = 2, 5 and 
10 GeV 2 . The curves are predictions from the statistical approach. 
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4 Inclusive neutral and charged current e p 
cross sections 



The neutral current DIS processes have been measured at HERA in a kine- 
matic region where both the 7 and the Z exchanges must be considered. The 
cross sections for neutral current can be written, at lowest order, as 



d 2 rr ± 
u °NC 

dxdQ 2 
where 



2na 
~x~Q A 



2 r 



Y + F 2 (x, Q 2 ) t Y_xF 3 (x, Q 2 ) - y 2 F L (x, Q 2 ) 



(17) 



F 2 (x,Q 2 



FT 1 - VeXz(Q 2 )G 2 (x, Q 2 ) + (a 2 + v 2 e ) X 2 z (Q 2 )H 2 (x, Q 2 ) 



;is) 



xF 3 (x, Q 2 ) = -a eXz (Q 2 )xG 3 (x, Q 2 ) + 2a e v eX 2 z (Q 2 )xH 3 (x, Q 2 ) . (19) 

The structure function Fl(x,Q 2 ) is sizeable only at high y and the other 
structure functions introduced above, have the following expressions in terms 
of the parton distributions 

[Fr, G 2 , H 2 ] (x, Q 2 ) = [Q)> 2 Qf v f> a ) + v f] Q 2 ) + Q 2 )) , 

f 

[xG 3 ,xh 3 ](x,q 2 ) = J2l 2 Qf a f> 2a f v f}( X( iA x 'Q 2 )- X( iA x >Q 2 ))- (2°) 

/ 



Here the kinematic variables are y = Q 2 /xs, Y± = 1±(1 — y) 2 , yfs = -J E e E p , 
E e and E p are the electron (positron) and proton beam energies respectively. 
Morever, Vi and are the vector and axial-vector weak coupling constants 
for the lepton e and the quark /, respectively, and Qf is the charge. The 
function Xz{Q 2 ) is given by 



4 sin 2 6 W cos 2 6 W Q 2 + M 2 



X,(Q 2 ) = 7—2-. — ^ ,,2 , (21) 



where 9w is the weak mixing angle and Mz is the Z-boson mass. The reduced 
cross sections are defined as 
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Our predictions are compared with HI and ZEUS data in Figs. Eland[TUJ as 
a function of x, in a broad range of Q 2 values and the agreement is excellent. 



The charged current DIS processes have been also measured accurately 
at HERA in an extented kinematic region. It has a serious impact on the 
determination of the unpolarized parton distributions by allowing a flavor 
separation because they involve only the W ± exchange. The cross sections 
are expressed, at lowest order, in terms of three structure functions as follows 



Al^cc (~<2 

a °Barn _ ^ F IV1 W 



dxdQ 2 4tt (Q 2 + M\ 



2 \2 

w) 



[Y + F~(x,Q 2 )-y 2 Fl c (x,Q 2 ) 



+Y_xF~(x,Q 2 )} , (23) 
and the reduced cross sections are defined as 



a cc (x,Q 2 ) 



L4vr {Q 2 + M 2 V ) 2 



~ l d 2 a cc 



dxdQ 2 

At leading order for e~p — > u e X with a longitudinally polarized beam 



(24) 



F?(x,Q 2 ) = x[u(x,Q 2 ) + c(x,Q 2 ) + d(x,Q 2 )+s(x,Q 2 )] 
xF^ c (x,Q 2 ) = x[u(x,Q 2 ) + c(x,Q 2 )- d(x,Q 2 ) - s(x,Q 2 )} , (25) 

and for e + p — > v e X 

F 2 cc (x,g 2 ) = x[d(x,Q 2 ) +s(x,Q 2 ) +u(x,Q 2 ) +c(x,Q 2 )} 
xFl c {x,Q 2 ) = x[d(x,Q 2 ) +s(x,Q 2 ) -u(x,Q 2 ) -c(x,Q 2 )} . (26) 

At NLO in QCD F£ c is non zero, but it gives negligible contribution, except 
at y values close to 1. Our predictions are compared with HI and ZEUS data 
in Figs. ITTl [T2l [T3l andfT4l as a function of x in a broad range of Q 2 values and 
vice versa. The agreement is very good, but unfortunately since the highest 
x value is only 0.42, it does not allow to clarify the situation regarding the 
large x behavior, as already noticed above. 
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Figure 9: The reduced neutral current cross section a, as a function of 
x, for different fixed values of Q 2 . Reaction e~p at s/s = 320GeV, e + p at 
\/s = 319GeV. Data from HI ^3123]. The curves are predictions from the 
statistical approach. 
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Figure 10: Same as Fig. El Data from ZEUS [201 - 
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Figure 12: The reduced charged current cross section a, in e p reactions as 
a function of Q 2 , for different fixed values of x. Data from HI [TBI I25j. 
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Figure 14: The reduced charged current cross section a, in e p reactions as 
a function of Q 2 , for different fixed values of x. Data from ZEUS \'27\. 
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5 Charged current neutrino cross sections 



The differential inclusive neutrino and antineutrino cross sections have the 
following standard expressions 



d 2 a^ 
dxdy 



G 2 F M p E u 

7T(1 + '" 



M, 



xy 2 F^\x, Q 2 ) + (1 - y - Mn^L )F ^), 



2E V 



y 



±(y-» Y )xF^\x,Q 2 ) 



x,Q 2 ) 
(27) 



y is the fraction of total leptonic energy transfered to the hadronic system 
and E u is the incident neutrino energy. F 2 and F 3 are given by Eq. (|25|) for 
vp and Eq. (|26|) for up, and F\ is related to F 2 by 



2xF x 



1 + 4M n 2 i 2 

—F 2 

1 + R 



(28) 



where R = ctl/ct, the ratio of the longitudinal to transverse cross sec- 
tions of the W-boson production. The calculations are done with sin 2 9w = 
0.2277 ±0.0013 ±0.009 obtained by NuTeV [2HI and the comparison with the 
CCFR and NuTeV data are shown in Fig. ED As expected, for fixed x, the 
y dependence is rather flat for neutrino and has the characteristic (1 — y) 2 
behavior for antineutrino. 



This can be extrapolated to evaluate the cross section of ultrahigh energy 
neutrinos with nucleons. The total cross section at a given neutrino energy 
reads 

Our prediction for this total charged current cross section, for an isoscalar nu- 
cleon N = l/2(p±n), versus the neutrino energy, is displayed in Fig. Upland 
it has the expected strong energy increase. We have not calculated the cor- 
responding neutral current cross section, which is known to be a factor three 
or so smaller. This new information is certainly valuable to the large scale 
neutrino telescopes, for the detection of extraterrestrial neutrino sources. 
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Figure 15: Differential cross section v{y)N for E v = 85GeV, as a function of 
y. Data are from CCFR [2E] (white circles) and NuTeV experiments [HP] EH] 
(black circles). 
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6 Drell-Yan dilepton cross sections 



A very important source of information for q(x) distributions comes from 
Drell-Yan dilepton processes, whose cross sections are proportional to a com- 
bination of products of q(x) and q(x) distributions. The cross section ody{j>P) 
for pp — > at the lowest order, has the simplified form 



M- 



d 2 a DY (pp) 
dMdxp 



lira 



9(xi + x 2 



e i [toMtoM + QiMQiM] 



(30) 



where M is the invariant mass of the produced muon pair, x\ and x 2 refer to 
the beam and target respectively, xp = X\ — x 2 and M 2 = Xix 2 s, where \fs 
is the center of mass energy of the collision. Clearly at NLO one should add 
the Compton processes contributions to the above qq annihilation terms. 
More recently the NuSea Collaboration has released the data on the Drell-Yan 
cross sections <Tdy(pp) and aoY^pd) for proton-proton and proton-deuterium 
collisions at 800 GeV/c They are displayed in Fig. [T7| as a function 

of xp for selected M bins, together with our predictions. The agreement is 
fairly good, mainly in the small mass region, but in order to evaluate it more 
precisely, we have plot in Fig.^Jthe ratios of experiment versus theory, using 
a broader set of data. 

Let us now come back to the extraction of the ratio d/u from these data. For 
large xp, namely X\ » x 2 and small M, we have 



2a DY (pp) 



1/2 



d(x 2 ) 
u(x 2 ) 



(31) 



so the measurement of this cross sections ratio is directly related to d(x)/u(x) 
for small x. For large x one needs to use small xp and large M values and 
we have now for x\ ~ x 2 



2a DY {pp) 



1/2 



, ^d(x)_ , ^d(x)_ , 2 d ( x ) d ( x ) 
u(x) u(x) u(x) u(x) 

u(x) u(x) 



(32) 



Therefore the falloff at large x oi ar>Y(pd)/2ap,Y(pp) observed in Ref.fSU] can- 
not be directly related to the falloff of d(x)/u(x), since d(x)/u(x) is also de- 
creasing for large x, as shown previously (see Fig.0J). The use of Eq. (}3*T]) will 
lead to an underestimation of d(x)/u(x). We also notice in Fig. Elan exper- 
imental point for (Jdy{pp) in the bin with M in the range (10.85,11.85)GeV 
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Figure 17: Drell-Yan cross sections per nucleoli at ^/s = 38.8GeV for pp and 
pd as a function of xp for selected M bins. Solid curve pp, dashed curve pd. 
Experimental data from FNAL E866 [32] • 

at xf — 0.05, two standard deviations above our curve, which might very 
well be one of the reason for the dramatic falling off of d(x)/u(x) for x ~ 0.3, 
reported by NuSea. Obviously more accurate data in this region are badly 
needed. 
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7 Single jet and n° inclusive productions 

A precise determination of parton distributions allows us to use them as 
input information to predict strong interaction processes, for additional tests 
of pertubative QCD and also for the search of new physics. Here we shall 
test our statistical parton distributions for the description of two inclusive 
reactions, single jet and 7r° productions. The cross section for the production 
of a single jet of rapidity y and transverse momentum p?, in a pp collision is 
given by 



dp 3 ~f 1 + $ij ft Jx 2x a — xre y 

_ d/T,.. 1 

(33) 



at 



where xt = 2pt/\/s, x = xxe y /(2 — XTe~ y ), Xb = x a xxe~ y / \2x a — xxe y ) 
and y/s is the center of mass energy of the collision. In the above sum, 
i,j stand for initial gluon-gluon, quark-gluon and quark-quark scatterings, 
d&ij/dt are the corresponding partonic cross sections and Q 2 is the scaling 
variable. The NLO QCD calculations at C(af), were done using a code 
described in Ref. jHH], based on a semi-analytical method within the "small- 



cone approximation" 7 . 

In Fig. El our results are compared with the data from CDF and DO 
experiments (HHl EZI- Our prediction agrees very well with the data up to 
the highest Ex (or p?) value and this is remarkable given the fact that the 
experimental results are falling off over more than six orders of magnitude, 
leaving no room for new physics. For completeness, we also show in Fig. 1201 
the DO data, for several rapidity bins, using a (Data-Theory)/Theory pre- 
sentation. 

Next we consider the cross section for the inclusive production of a 7r° 
of rapidity y and transverse momentum pr, in a pp collision, which has the 

7 We thank Werner Vogelsang for providing us with the numerical values, resulting from 
the use of our parton distributions 
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following expression 



E n d 3 a/dpl = j dx a dx b f a/p (x a ,Q 2 ) x 

abc 

fb/p(x b , Q^ ArO/c^c, Q 2 ) d&/d ^ ab ^ cX) ; (34) 
7TZ C 

where the sum is over all the contributing partonic channels ab — *> cX and 
da/dt is the associated partonic cross section. In these calculations f a / p , fb/ p 
are our parton distributions and D n o/ c is the pion fragmentation function. 
Our calculations are done up to the NLO corrections, using the numerical 
code INCNLL of Ref. [HH] and for two different choices of fragmentation 
functions namely, BKK of Ref. |1U] and KKP of Ref. |41j . and we have 
checked that they give similar numerical results. We have compared our 
predictions to two different data sets at ^Js = 200GeV from PHENIX and 
STAR at RHIC-BNL. The results are shown in Figs. |2U and [221 and the 
agreement is good, both in central rapidity region (PHENIX) and in the 
forward region (STAR). This energy is high enough to expect NLO QCD 
calculations to be valid in a large rapidity region, which is not the case for 
lower energies |4*lj . 
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Figure 20: Comparison between the statistical model and the DO [SB] single 
jet cross sections in pp at y/s = 1.8TeV, as a function of E T and rapidity r\. 
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Figure 21: Inclusive 7r° production in pp reaction at y/s = 200GeV as a func- 
tion of Pt, scale \i = p?- Data from PHENIX [32] ■ Solid curve fragmentation 
functions from KKP |41j . dashed curve from BKP |4()j . 
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Figure 22: Inclusive 7r° production in pp reaction at yfs = 200GeV as a 
function of E v . Data from STAR [43]. Solid curve fragmentation functions 
from KKP jlT] , dashed curve from BKP [10] . 
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8 Concluding remarks 

We have shown that this simple approach of the statistical parton distri- 
butions provides a good description of recent data on unpolarized and po- 
larized DIS and on several hadronic processes. Since it involves only eight 
free parameters, we have tried to relate them to some specific properties 
of the parton distributions, but we do not have yet a full understanding of 
their physical interpretation. It is important to stress that we have simul- 
taneously the unpolarized and the polarized parton distributions, which is a 
unique situation. The main features of our distributions agree with other sets 
available in the literature, both in the unpolarized case |H3 HB] see Fig. 1221 
and in the polarized case |T71 HH1 HI] ■ We show in Fig. 1221 a comparison with 
MRST and CTEQ, where one observes that the essential differences lie in 
the small x region. We have also identified some physical observables and 
kinematic regions, where we can make definite predictions. In particular, let 
us recall a slow decreasing behavior of d(x)/u(x) for x > 0.6, the fact that 
d(x)/u(x) should remain larger than one for x > 0.3, the signs Au(x) > 
and Ad(x) < and our choice for AG(x). All these are real challenges and 
we look forward to new precise experimental data in the future. 
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